Pairing and condensation in a resonant Bose-Fermi mixture 
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We study by diagrammatic means a Bose-Fermi mixture, with boson-fermion coupling tuned by a 
Fano-Feshbach resonance. For increasing coupling, the growing boson-fermion pairing correlations 
progressively reduce the boson condensation temperature and make it eventually vanish at a critical 
coupling. Such quantum critical point depends very weakly on the population imbalance and for 
vanishing boson densities coincides with that found for the polaron-molecule transition in a strongly 
unbalanced Fermi gas, thus bridging two quite distinct physical systems. 
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One of the main reasons behind the recent great in- 
terest in ultracold gases, is the possibility to reproduce 
physical systems relevant to other areas in physics, with 
a flexibility and a degree of tunability of physical parame- 
ters which is unimaginable in the original system of inter- 
est. The very same flexibility can also be used, however, 
to construct novel physical systems. A noticeable exam- 
ple of such novel systems, namely, resonant Bose-Fermi 
mixtures, will be here at issue. We will be interested 
in particular in the competition between boson-fermion 
pairing correlations and bosonic condensation, occuring 
in a mixture of single-component bosons and fermions 
when the boson-fermion pairing is made progressively 
stronger by means of a Fano-Feshbach resonance. 

Initial theoretical studies of ultracold Bose-Fermi mix- 
tures considered mainly non- resonant systems [H-Q) 
where boson-fermion pairing is irrelevant, and studied, 
within mean-field like treatments, the tendency towards 
collapse or phase separation (as motivated by the first 
experimental results on non-resonant Bose-Fermi mix- 
tures @). Bose-Fermi mixtures in the presence of a Fano- 
Feshbach resonance have been subsequently considered 
in mostly for narrow resonances 0, 1(3- HI EH 

and/or in specific contexts such as optical lattices, re- 
duced dimensionality, zero temperature or vanishing den- 
sitiy of one component [1, Ei > HI 17 1 . 

Most of current experiments on Bose-Fermi mix- 
tures 19|-23| appear however to be closer to the case of 



body calculations (beyond mean- field). 

We thus consider the following (grand-canonical) 
Hamiltonian: 



a broad resonance, which is characterized by the small- 
ness of the effective range parameter ro of the interac- 
tion potential with respect to both the average inter- 
particle distance and the scattering length [18[. Un- 
der these conditions, the resonant Bose-Fermi mixture 
is accurately described by a minimal Hamiltonian, made 
just by bosons and fermions mutually interacting via an 
attractive point-contact potential. This Hamiltonian is 
simpler and more "fundamental" in character than its 
counterpart for a narrow resonance because of the ab- 
sence of any other parameter besides the strength of 
the boson-fermion interaction. In addition, its simplicity 
allows the implementation of more sophisticated many- 
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Here ipt ( r ) > creates a particle of mass m s and chemical 
potential fi a at spatial position r, where s = B,F indi- 
cates the boson and fermion atomic species, respectively, 
while wo is the bare strength of the contact interaction 
(we set h = fcs = 1 throughout this paper). Ultraviolet 
divergences associated with the use of a contact interac- 
tion in ([T]) are eliminated, as for two-component Fermi 
gases [24j , by expressing the bare interaction vq in terms 
of the boson-fermion scattering length a via the (diver- 
gent) expression 1/vq = m r /(2na) — /2m r /k 2 (ik/(27r) 3 , 
where m r = msmp / '{tub + mp) is the reduced mass. 

The Hamiltonian (JTJ does not contain explicitly the 
boson-boson interaction. Provided it is repulsive and 
non- resonant, this interaction yields in fact a mean- field 
shift, which consists in a simple redefinition of the bo- 
son chemical potential (attractive boson-boson interac- 
tions would lead to mechanical instability and are ex- 
cluded from our consideration). S-wave interaction be- 
tween fermions is finally excluded by Pauli principle. 

The effective coupling strength in the many-body sys- 
tem is determined by comparing the boson-fermion scat- 
tering length a with the average interparticle distance 
n^ 1 / 3 (with n = ub + np, where and np are the bo- 
son and fermion particle number density, respectively). 
In particular, we will use the same dimensionless cou- 
pling strength (kpa)" 1 normally used for two-component 
Fermi gases, where the wave- vector kp = (S^n) 1 ^ (note 
that kp coincides with the noninteracting Fermi wave- 
vector k F = {%'K 2 nF) 1 / z only for = np). 

The expected behavior of the system is clear in the 
two opposite limits of the boson-fermion coupling. In 
the weak-coupling limit, where the scattering length a 
is small and negative (such that (fcpa) -1 <C —1), the 
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FIG. 1: T-matrix diagrams for the fermionic and bosonic self- 
energies in the normal phase. Full lines represent bare bosonic 
(BB) and fermionic (FF) Green's functions. Broken lines rep- 
resent bare boson- fermion interactions Vq. 



two components will behave essentially as ideal Bose and 
Fermi gases. Bosons will condense at T c = 3.31n^ 3 /tub, 
while fermions will fill out (at sufficienty low tempera- 
ture) a Fermi sphere with radius k F , the chemical po- 
tentials hb and \i F being modified with respect to their 
noninteracting values by the mean-field shifts 2irn F a/m r , 
and 2irriBa/m r , respectively. In the opposite strong- 
coupling limit where a is small and positive (such that 
(fcpa) -1 ^> 1 ), bosons will pair with fermions to form 
tightly bound fermionic molecules, with binding energy 
eo = l/(2m r a 2 ). In particular, in systems where tib < 
rip, on which we focus in the present paper, all bosons 
will eventually pair with fermions into molecules, sup- 
pressing condensation completely. 

In a complete analogy with the BCS-BEC crossover 
problem in a two-component Fermi gas (2~3 |. the choice 
of the self-energy diagrams will be guided by the criterion 
that a single set of diagrams should recover the correct 
physical description of the two above opposite limits. In 
the normal phase above the condensation critical tem- 
perature T c , to which we restrict in the present paper, 
this condition is met (as we shall see below) by the T- 
matrix choice of diagrams represented in Fig. 1, yielding 
the expressions: 

dP 
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(2tt) 3 ^ 

for the bosonic and fermionic self-energies Y>b and Hp, 
where G° B and G° F are bare bosonic and fermionic Green's 
functions and T(P) is the many-body T-matrix: 

r(P) = _ r mr , r d P 



m r 

2^a + J (2tt) 3 

i-/Kf(p-p)] + &Kb(p)] 

£>(P-p)+fn (p)-fflm 

In the above expressions q = (q, w„), 
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FIG. 2: Critical temperature (in units of E F — k%/2m F ) 
for condensation of bosons vs the boson-fermion coupling 
(fcFts) -1 for different values of the density imbalance (n F — 
tib)/ n in a mixture with tub = m F . 



P = (P, f2 m ), where u) v = 2ttv/ /3 and ui n = (2n + l)n//3, 
m = (2m + l)7r/ /3 are bosonic and fermionic Matsub- 
ara frequencies, respectively (fi — 1/T being the inverse 
temperature and v, n, m integer), while f(x) and b(x) are 
the Fermi and Bose distribution functions and £ s (p) = 
p 2 / (2m s ) — n s . The above self-energies determine finally 
the dressed Green's functions Gj 1 = G° _1 — T, s entering 
the particle numbers equations: 

= E G -( k '-«)e- 0+ . (6) 

Equations (J2HU) fully determine the thermodynamic 
properties of the Bose-Fermi mixture in the normal 
phase. This phase is characterized by the condition 
Hb — £_b(<7 = 0) < 0; upon lowering the temperature, 
condensation will then occur when /is = Sb (q = 0). 

Analytic results can be obtained in the two opposite 
limits of the boson-fermion coupling [2^]. For weak- 
coupling T(P) ~ — ^p, yielding ~ 2im F a/m r and 
Y, F ~ 2TrriBa/m r , in accordance with the expected mean- 
field shifts. For strong-coupling T(P) gets instead pro- 
portional to the molecular propagator: 



r(p) 



2tt 



1 



mj.a m m - 



(7) 



(k,w„ 



where /jm = fJ-B+fJ-F+^o is the molecular chemical poten- 
tial and M = raj + m F . Insertion of this propagator in 
Eqs. (|2Hni) leads, for i%b < n F , to [iu — (67r 2 ns ) 2 / 3 /2M, 
Hf ~ [67r 2 (n F - n B )] 2 ^ 3 /2m F and fi B = Mm - - £o — 
— eo, as expected when all bosons pair with fermions. 

For intermediate values of the boson-fermion coupling 
Eqs. ([2]— [6J) need to be solved numerically. The results for 
the condensation critical temperature T c vs the boson- 
fermion dimensionless coupling (kpa) -1 are presented in 
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FIG. 3: Boson chemical potential at the critical temperature 
T c (in units of Ef) vs the boson-fermion coupling {kpa)~ 
for different values of the density imbalance (tif — riB)/n in a 
mixture with tub = rap. The corresponding fermion chemical 
potential is reported in the inset. 



FIG. 4: Boson (full line) and fermion (dashed-dotted) mo- 
menta distribution curves for ub = np, ms = m-F, 
{kpay 1 = 1.63 and T = T c ~ 0M5E F . See text for the 
meaning of the dashed curve. The inset reports the fraction 
of molecules vs (fe^ci) -1 for the same mixture at T = O.IEf- 



Fig. 2 for a mixture with tub = mp at several values of 
density imbalance. [26j j All curves start in weak-coupling 
from the corresponding noninteracting values, and de- 
crease monotonically for increasing coupling due to the 
growing pairing correlations, which tend to deplete the 
zero momentum mode and distribute the bosons over a 
vast momentum region, as required to build the internal 
molecular wave-function. The critical temperature van- 
ishes eventually at a critical coupling, corresponding to 
a quantum critical point which separates a phase with 
a condensate from a phase where molecular correlations 
are so strong to deplete the condensate completely. 

Remarkably, the critical coupling value depends very 
weakly on the degree of density imbalance: all curves 
terminate in the narrow region 1.6 < (kpa) -1 < 1.8. 
In this respect, it is interesting to consider the limit 
Tib 0, where the critical coupling can be calculated 
independently by solving the problem of a single boson 
immersed in a Fermi sea. This is actually the same as 
a spin-down fermion immersed in a Fermi sea of spin-up 
fermions, since for a single particle the type of statistics 
is immaterial. The critical coupling reduces then to that 
for the polaron-to-molecule transition, recently studied in 
the context of strongly imbalanced two- component Fermi 
gases 27-33J. In particular, the solution of our equations 
in this limit yields the critical coupling (kpa)^ 1 = 1.60, 
in full agreement with the value (kta)^ 1 = 1.27 (where 
k F = 2 1 / 3 kp) reported in [32], l33j ] for the polaron-to- 
molecule transition, when treated at the same level of ap- 
proximation 34|, |36[ . The meeting of the properties of a 
Bose-Fermi mixture with those of a two-component Fermi 
mixture at this point of the phase diagram is quite re- 
markable, especially because, according to our results of 
Fig. 1, this "universal" point sets the scale for the quan- 
tum phase transition for all boson densities ub < rip. 

The chemical potentials /xg and jip at T c vs. (kpa)~ l 



are reported in Fig. 3, for different values of the density 
imbalance. The two chemical potentials behave quite dif- 
ferently. The fermion chemical potential remains almost 
constant in the whole range of coupling considered; the 
boson chemical potential, on the other hand, diminishes 
quite rapidly with increasing coupling while depending 
little on the density imbalance. This different behav- 
ior results from the concurrence of several factors. For 
weak-coupling, the increasing (negative) mean-field shift 
of the fermion chemical potential for increasing coupling 
is partially compensated by the decrease of the tempera- 
ture when moving along the critical line. On the molecu- 
lar side, the fermion chemical potential is instead deter- 
mined by the Fermi energy of the unpaired fermions plus 
a mean-field shift caused by interaction with molecules. 



Pauli repulsion makes this interaction repulsive [8j, |25 | 
thus keeping the fermion chemical potential positive. The 
boson chemical potential, on the other hand, interpolates 
between the mean-field value 2impa/m r in weak cou- 
pling and /is ~ — eo m strong coupling, as required by 
molecule formation. 

Figure 4 reports the momentum distributions rts(k) 
and rii?(k) (as obtained before momentum integration 
in Eqs. ©-([H])) at T c , for a mixture with ub = Tip, 
at the coupling value {kpa)^ 1 = 1.63 (approximately 
at the quantum critical point). The two distributions 
are markedly different at low momenta, consistently with 
their different statistics, but coalesce into the same be- 
havior just after the step in the fermion momentum 
distribution. This common behavior corresponds to 
the function 7iM|0(k)| 2 (dashed line in Fig. 4), where 
<^>(k) = (87ra 3 ) 1 / 2 /(k 2 a 2 + 1) is the normalized two-body 
internal wave function of the molecules, while the coeffi- 
cient um represents their density. In the present case 
njvf = 0.89 ns, showing that a fraction of bosons re- 
mains unpaired but still does not condense even at such 
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a low temperature. The extrapolation of these results 
at exactly zero temperature indicates the existence (in 
a coupling range starting right after the quantum crit- 
ical point) of quite an unconventional Bosc liquid, cor- 
responding to the unpaired bosons, which do not con- 
dense even at zero temperature. The fraction of unpaired 
fermions is instead more conventional and consists in a 
Fermi liquid, which is responsible for the jump in the 
fermion momentum distribution. In particular, the posi- 
tion of the jump in Fig. 4 at |k| ~ 0.47 kp, corresponds 
to an "enclosed" density of 0.10 np, in good numerical 
agreement with the value 0.11 nj? obtained from the Lut- 
tinger theorem [37j ] for the fraction (np — Um) of un- 
paired fermions (using the value Um = 0.89 ns extracted 
independently above). Note finally that the number of 
unpaired bosons progressively decreases by increasing the 
coupling, as expected on physical grounds, reaching even- 
tually a 100% conversion of bosons into molecules, as 
shown in the inset of Fig. 4, where the ratio n^i/nB is 
reported vs coupling at a constant temperature. 

We wish to conclude finally by commenting on three- 
body losses and mechanical instabilities, that could pre- 
vent the experimental observation of the many-body 
physics described above. The choice of restricting our 
analysis to mixtures with a majority of fermions was 
precisely aimed at reducing the importance of these nui- 
sances. Three-body losses are dominated by processes 
involving two bosons and one fermion, with a rate propor- 
tional to n 2 B riF [2^| . They can thus be controlled by keep- 
ing riB sufficiently small with respect to np. A predomi- 
nance of fermions should also reduce the tendency to col- 
lapse, because of the stabilizing effect of Fermi pressure. 
As a matter of fact, we have verified that the compress- 
ibility matrix dfj, s /dn s > remained positive definite for all 
couplings, temperatures, and densities np > ns consid- 
ered in our calculations. Apparently, pairing correlations 
act to protect the system from the mean field instabili- 
ties dominating the phase diagram of non-resonant Bosc- 
Fermi mixtures. Resonant Bose-Fermi mixtures with a 
majority of fermions appear thus promising systems for 
the observation of a rich many-body physics. 
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